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FOREWORD 


This report was prepared under Work Unit 240410A1 entitled "Aerodynamic 
Design and Analysis Methods" by Howard T. Emsley of the Aerodynamic Methods 
Group, Aerodynamics and Airframe Branch, Aeromechanics Division, Flight 
Dynamics Directorate, Wright Laboratory, Wright-Patterson Air Force Base, 
Ohio. 

The report documents PLUTO (Poisson Laplace U TFI Orthogonal), a three 
dimensional gridding and smoothing program, which was developed in-house 
between July 1987 and May 1991. The coding was primarily performed by Capt. 
David Amdahl and Mr. William Strang of the Aerodynamic Methods Group. 

This document has been reviewed and approved. 





Howard T. Emsley 
Aerospace Engineer 


Dennis Sedlock 

Chief, Aerodynamics and Airframe Branch 
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INTRODUCTION 


D 


j 


) PLUTO (Poisson Laplace U TFl Orthogonal), is a three-dimensional gridding 
and smoothing program developed in-house at WL/FIMM^ Its purpose is to: (1) 
generate three-dimensional grids from the six two/dimensional faces of a 
computational grid, (2) check grid quality, (3) smooth grids to improve 
spacing, smoothness and orthogonality, and (4) output grids in user selected 
ascii or binary formats.' 




Initial grids are established by PLUTO with a three-dimensional algebraic 
transfinite interpolation (TFI) scheme. For many applications, where simple 
block shapes and spacings are used, an algebraic grid may be sufficiently 
smooth and orthogonal for use. As a check, QBERT (Quantitative Block ERror 
Tester) is run on each block to generate diagnostics about the grid’s 
spacing, skewness and right-handedness. —- 


If the user is unsatisfied with the initial grid, Laplace and/or Poisson 
smoothing techniques are available, and orthogonality constraints can be 
enforced separately on each block face. By handling each grid block 
independently the user is afforded great flexibility and control over the 
finished grid. For example, complex block shapes whose initial algebraic 
grid are often crossed, can be smoothed with the Laplace technique to 
uncross the grid. 


Once a finished grid block is generated, PLUTO also assists in outputting 
the blocks in a final format. This includes the capabilities to cut large 
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blocks into smaller blocks (this is used to subdivide large computational 
blocks to run on machines with limited memory), and/or output the grid in 
one of nine different formats. 


BACKGROUND 

Development of PLUTO has been somewhat of an evolutionary process. In 1987, 
the Aerodynamic Methods Group (WL/FIMMA) began a concerted effort to develop 
rapid geometric modeling and grid generation capabilities for the Air Force. 
Programs such as I3G/VIRG0^, PLUTO and GRIDGEN^ were the results of in-house 
and contracted efforts between 1987 and 1991. Each code has its specific 
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applications where it excels. I3G was originally developed under contract 
as surface modeller and then it was modified in-house (VIRGO) to support 
grid generation needs. PLUTO was then developed to generate three- 
dimensional grids as discussed above. Concurrently, GRIDGEN (GRIDBLOCK, 
GRIDGEN2D,and GRIDGEN3D) was developed under contract to generate grids 
around predefined surface databases. GRIDGEN3D performs the same function 
in the GRIDGEN grid generation process as PLUTO does for I3G/VIRG0. 

PLUTO is not a revolutionary piece of coding, but its useability, 

flexibility, and inclusion of diagnostics (QBERT) are the things which make 

this code noteworthy ana in high use. PLUTO incorporates the work of a 

variety of people in the grid generation field. The TFI method used is 

attributed to Dr. Soni, the Poisson control functions are from Dr. Joe 
5 

Thompson’s work , and the orthogonal control functions are a variation of 

6 

those found in Sorenson’s GRAPE Method . 

PLUTO is written primarily in standard FORTRAN-77 however a few routines 
have been written in C. The code was originally written to operate on a 
Cray X-MP, however the version documented here was tailored for operation on 
an Iris workstation. Portability to other machines is expected to be 
straight forward, however, memory constraints and machine speed should be 
considered when porting to a less advanced host. 
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PROGRAM EXECUTION 


When PLUTO is executed, the user is prompted to provide four different 
files: Input Data File, Output Data File, Input Grid File and the Output 
Grid File. The input files are existing files that the user will provide 
and the output files are new files the program will create. 

The Input Data File is the equivalent of a job file on a Cray. It 
provides the program with all the flags necessary to direct the needed 
operations. Everything from block dimensions to convergence residuals 
are specified in this file. (See examples) 

The Output Data File is the file where QBERT will write diagnostics 
about the blocks in the grid. The user should examine this file to 
determine the results of the run. At a minimum, this file will contain 
each block’s dimensions and information on the worst cells of the 
block. Additional information is provided when required 

The Input Grid File is the file on which PLUTO will operate. This file 
can be in ascii or binary format and will provide either the six 
starting faces of each block or a three-dimensional grid for each 
block. (See examples for further details) 

The Output Grid File will contain the new three-dimensional grid. The 
format of this file is controlled with flags located in the Input Data 
File. 


INPUT DATA FILE 

PLUTO can be used to perform a wide variety of different functions for 
the user. Control of the program is handled almost exclusively in the 
Input Data File. A general explanation of this file follows, but often 
the examples included may provide the best guide. 


3 





Lines 1 and 2 of the file consist of a comment line explaining the 
inputs, a value for smoothing relaxation (values between 0. and 2.), 
the maximum number of iterations the smoother will run, and how often 
the screen will display the convergence information (see below). The 
solution can be over or under relaxed by varying the relaxation range. 

♦♦RELAXATION NO. OF ITERATIONS PRINT OUT/ITERATIONS (line 1) 

.8 50 5 (line 2) 

Lines 3 and 4 consist of a comment line labeling the input, a flag for 
the convergence (average change of the block or maximum change of any 
one point), and the residual value for a converged solution. 

♦ ♦CONVERGENCE(0=AVER, 1=MAX. ) FINAL RESIDUAL (line 3) 

1 .000001 (line 4) 

Lines 5 thru 10 contain four comment lines and two lines containing 

form and format flags for the grid files. File form refers to ascii or 

binary format, and file format refers to the geometric format of the 
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input and output files. PANAIR and BLOCK formats provide only the six 
faces of each block while GRIDGEN, MERCURY 8 , TEAM 9 and PL0T3D 10 expect 
a full three-dimensional grid. The user is directed to each code’s 
manuals for format specifics. 


♦♦FORM (1=BINARY,2=ASCII) (line 5) 

♦ ♦FORMATS (1=PANAIR,2=BL0CK,3=GRIDGBN,4-MERCURY,5=TEAM,6=PL0T3D) (line 6) 
♦♦INPUT FILE FORM INPUT FILE FORMAT (1,2, AND 3 ARE ASCII ONLY)(line 7) 

2 1 (line 8) 

♦ ♦OUTPUT FILE FORM OUTPUT FILE FORMAT (ONLY 4,5, AND 6) (line 9) 

1 6 (line 10) 


Lines 11 thru 16 contain information for the orthogonality conditions. 
The columns correspond to P, Q, and R while the rows correspond to 
sides 1 through 6. The user need not modify these values, unless the 
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solution will not converge. If this occurs, the user is instructed to 
raise the values of all P, Q and R values and lower the relaxation 
value on line 2. The values shown below are used as exponential 
exponents, therefore values less than or equal to zero should not be 
used. (Enforcing orthogonality in three dimensions is a complicated 
task, and obtaining a converged solution may require a healthy dose of 
experimentation. In general, there is no standard "fix" to obtain 
convergence.) 


.65 

.25 

.25 

(line 11) 

.65 

.25 

.25 

(line 12) 

.25 

.65 

.25 

(line 13) 

.25 

.65 

.25 

(line 14) 

.25 

.25 

.65 

(line 15) 

.25 

.25 

.65 

(line 16) 


Lines 17 through 20 contain three comment lines and the number of 
blocks the program should be reading from the Input Grid File. 

Starting on line 21, the Input Data File should contain information 
about the blocks in question. A variety of different inputs can appear 
after line 20 depending on the application desired. An unrealistic 
combination of options is shown below to illustrate the kinds of 
information that may be required for each block. 

In general, PLUTO will require the block dimensions, the mode and a set 
of orthogonality flags for each block. The mode selected will 
determine the operation that will be performed on the block, and the 
orthogonality flags will determine which sides will have orthogonality 
imposed (0 = no orthogonality, 1 = orthogonality). These flags will 
only be taken into account under modes 2, 3 and 4. Additional 
information will be required for blocks the user wishes to cut into two 
or more smaller blocks. 

When the grid is to be used on machines with limited memory, block 
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cutting is important to subdivide large computational grid blocks into 
smaller blocks. To cut a block, a ten should be added to the mode 
selected (e.g. mode 2 becomes mode 12). Following the mode line, the 
Input Data File should provide the number of new blocks and lines 
containing the looping indices of each new block (see block 3 below). 

If a block is being smoothed and cut, the block will not be cut until 
after the smoothing has been completed. By performing the functions in 
this order, slope continuity of the grid will occur between the newly 
cut blocks. 

The first block (lines 21, 22 and 23) has dimensions of 94 103 35 and 
will be generated with a Transfinite Interpolation scheme (mode 1). 

This format is generally used to generate the initial grid from the six 
block faces. The orthogonality flags in this case are of no 
consequence. 

The second block (lines 24, 25 and 26) with dimensions of 65 60 27 will 
be read in and written out with no manipulation of the grid. This 
assumes that the grid has already been generated as done in the 
previous block or possibly by another 3-D grid generator (such as 
GRIDGEN3D). This option can also be used if the user wants to rewrite 
existing blocks in another form or format without smoothing (see lines 
5 thru 10 for further details). Again the orthogonality flags are of 
no consequence. 

The third block (lines 27 to 29) with dimensions of 33 44 27 will be 
iterated on by a Poisson smoother for the number of iterations 
specified in line 2. Line 29 contains the six flags that can be set to 
impose orthogonality conditions at any or all of the block’s six faces 
(0 = no orthogonality, 1 = orthogonality). 

The fourth block (lines 30 to 35) with dimensions of 33 4 27 will be 
iterated on by a Laplacian smoother with orthogonality imposed on faces 
3, 5 and 6. After the defined iterations (line 2) or convergence, the 
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block will be cut into two smaller blocks at i=10. In this case, line 
30 contains the current dimensions of the block, line 31 indicates the 
smoother mode (2) and that a cut will take place (+10), line 32 
indicates number of new blocks, lines 33 and 34 provide the new 
dimensions, and line 35 provides the orthogonality flags (used before 
cut) . 


♦♦NUMBER OF BLOCKS 
**DDffiNSION 

♦ ♦MODE (0=READ/WRITE, 1=TFI, 2=LAPL, 3=PQIS, 4=P0IS w/CURV) 
4 

94 103 35 
1 

0 0 0 0 0 0 
65 60 27 
0 

0 0 0 0 0 0 
33 44 27 

3 

1 1 0 0 0 0 
33 4 27 

12 
2 

1 10 1 4 1 27 

10 33 1 4 1 27 

0 0 10 11 


(line 

17) 

(line 

18) 

(line 

19) 

(line 

20) 

(line 

21) 

(line 

22) 

(line 

23) 

(line 

24) 

(line 

25) 

(line 

26) 

(line 

27) 

(line 

28) 

(line 

29) 

(line 

30) 

(line 

31) 

(line 

32) 

(line 

33) 

(line 

34) 

(line 

35) 


7 





MODES 


PLUTO provides the user with 5 modes of operation: 

0=READ/WRITE 
1=TFI 
2=LAPLACE 
3=P0ISS0N 

4=P0ISS0N with CURVATURE 

The simplest of these is READ/WRITE which is appropriately named. This 
mode will cause the program to read a block in and write it back out 
without modifying the original block. It can be used quite effectively 
to transform one file format or form into another. 

TFI (Transfinite Interpolation) is an algebraic grid generator, which 
is often used to generate the original grid from the six block faces. 
For many applications, algebraic grids will be sufficiently smooth and 
no smoothing will be required. 

The Laplace mode tends to equalize the volume of the cells in a block. 
It can also be used to uncross algebraic grids that may be crossed due 
the block shape. Crossed grids must be uncrossed before the other 
smoothing modes are used. 

The Poisson and Poisson with curvature modes are also used to 
elliptically smooth the grid. In these modes, Thomas-Middlecoff 
control functions (based on the boundary arclengths) are linearly 
interpolated across the grid to control the spacing. In the latter 
mode, the control functions are modified by the curvatures of the 
adjacent boundaries. 

Only the three smoothers (modes 2,3 and 4) allow the user to control 
orthogonality at the block faces. 
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Experimentation and practice will provide the user a much clearer 
understanding of each modes applications. Two dimensional forms of 
these smoothers are used in I3G/VIRG0, and the user is directed there 
to experiment and view the different results that can occur. 
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EXAMPLES 


The following three example have been included to illustrate a variety of 
different Input Data Files. 

EXAMPLE 1 


In this example, we will generate four blocks of a grid from the six 
faces of each block in an ascii PANAIR formatted file. The three- 
dimensional grid will be output in a binary PL0T3D format and the grid 
will be generated algebraically with TFI. The Input Data File is as 
follows. 


★♦RELAXATION NO. OF ITERATIONS PRINT OUT/ITERATIONS 


.8 

50 

5 

*★CONVERGENCE(0=AVER, 

1=MAX. ) 

FINAL RESIDUAL 

1 


.000001 

★★FORM (1=BINARY,2=ASCII) 

★★FORMATS (1=PANAIR,2=BL0CK,3= 

GRIDGEN,4=MERCURY,5=TEAM,6=PL0T3D) 

★★INPUT FILE FORM 

INPUT FILE FORMAT (1,2, AND 3 ARE ASCII ONLY) 

2 

1 


★★OUTPUT FILE FORM 

OUTPUT FILE FORMAT (ONLY 4,5, AND 6) 

1 

6 


.65 

.25 

.25 

.65 

.25 

.25 

.25 

.65 

.25 

.25 

.65 

.25 

.25 

.25 

.65 

.25 

.25 

.65 


★★NUMBER OF BLOCKS 
★★DIMENSION 

★★MODE (0=READ/WRITE,1=TFI,2=LAPL,3=P0IS,4=PQIS w/CURV) 
4 
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94 103 35 

1 

0 0 0 0 0 0 
65 60 27 

1 

0 0 0 0 0 0 
33 44 27 

1 

0 0 0 0 0 0 
33 4 27 

1 

0 0 0 0 0 0 


EXAMPLE 2 

Now that a three-dimensional grid exists we will smooth blocks 2, 3 and 
4 for a maximum of 150 iterations while leaving block 1 as is. The 
Input Grid File is now a binary PL0T3D file, and the Grid Output File 
is an ascii MERCURY file. No orthogonality is specified. 


♦♦RELAXATION NO. OF ITERATIONS PRINT OUT/ITERATIONS 

.8 150 5 

*♦CONVERGENCE(0=AVER, 1=MAX. ) FINAL RESIDUAL 
1 .000001 
♦♦FORM (1=BINARY,2=ASCII) 


♦♦FORMATS (1=PANAIR,2=BL0CK,3=GRIDGEN,4=MERCURY,5=TEAM,6=PL0T3D) 
♦♦INPUT FILE FORM INPUT FILE FORMAT (1,2, AND 3 ARE ASCII ONLY) 
1 6 

♦♦OUTPUT FILE FORM OUTPUT FILE FORMAT (ONLY 4,5, AND 6) 

2 4 


.65 .25 .25 
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65 

.25 

.25 

25 

.65 

.25 

25 

.65 

.25 

25 

.25 

.65 

25 

.25 

.65 


**NUMBER OF BLOCKS 
♦♦DIMENSION 

♦♦MODE (0=READ/WRITE,1=TFI,2=LAPL,3=POIS,4=POIS w/CURV) 
4 

94 103 35 

0 

0 0 0 0 0 0 
65 60 27 

2 

0 0 0 0 0 0 
33 44 27 

3 

0 0 0 0 0 0 
33 4 27 

2 

0 0 0 0 0 0 


EXAMPLE 3 

We will now assume that the grid is fully smoothed. Block 1 however is 
too big for our solver therefore it must be cut into three pieces. The 
Input Grid File is now an ascii MERCURY file, and the Grid Output File 
will be an ascii MERCURY file. 


♦♦RELAXATION NO. OF ITERATIONS PRINT OUT/ITERATIONS 
.8 50 5 
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♦♦CONVERGENCE(0=AVER, 1=MAX. ) 


FINAL RESIDUAL 


1 .000001 
♦♦FORM (1=BINARY,2=ASCII) 

♦*FORMATS (1=PANAIR,2=BL0CK,3=GRIDGEN,4=MERCURY,5=TEAM,6=PL0T3D) 


♦♦INPUT FILE FORM 

INPUT 

FILE FORMAT (1,2, AND 3 

ARE ASCII ONLY) 

2 


4 



♦♦OUTPUT FILE FORM 

OUTPUT 

FILE FORMAT (ONLY 4,5, 

AND 6) 

2 


4 



.65 

.25 


.25 


.65 

.25 


.25 


.25 

.65 


.25 


.25 

.65 


.25 


.25 

.25 


.65 


.25 

.25 


.65 



♦♦NUMBER OF BLOCKS 
♦♦DIMENSION 

♦♦MODE (0=READ/WRITE,1=TFI,2=LAPL,3=P0IS,4=P0IS w/CURV) 


4 

94 103 35 

10 
3 


1 

94 

1 

30 

1 

35 

1 

94 

30 

60 

1 

35 

1 

94 

60 

103 

1 

35 

0 0 

0 0 

0 0 




65 

60 

27 





0 

0 0 0 0 0 0 
33 44 27 

0 

0 0 0 0 0 0 
33 4 27 

0 

0 0 0 0 0 0 
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